Laser cluster interaction in ambient magnetic fields for accelerating electrons in two stages without external injection

In the few-cycle pulse regime of laser-cluster interaction (intensity \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$>10^{16}\,\text{ W/cm}^{2}$$\end{document}>1016W/cm2, wavelength \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$> 780$$\end{document}>780 nm), laser absorption is mostly collisionless and may happen via anharmonic resonance (AHR) process in the overdense (cluster) plasma potential. Many experiments, theory and simulation show average absorbed energy per cluster-electron (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathcal {E}_A}$$\end{document}EA) close to the electron’s ponderomotive energy (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$U_\mathrm {p}$$\end{document}Up) in the collisionless regime. In this work, by simple rigid sphere model (RSM) and detailed particle-in-cell (PIC) simulation, we show enhanced \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathcal {E}_A}\approx$$\end{document}EA≈ 30–70\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$U_\mathrm {p}$$\end{document}Up—a 15–30 fold increase—with an external (crossed) magnetic field near the electron-cyclotron resonance (ECR). Due to relativistic mass increase, electrons quickly deviate from the standard (non-relativistic) ECR, but time-dependent relativistic-ECR (RECR) happens which also contributes to enhanced \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathcal {E}_A}$$\end{document}EA. Here laser is coupled to electrons in two stages, i.e, AHR and ECR/RECR. To probe further we retrieve the phase-difference \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Delta \psi$$\end{document}Δψ between the driving electric field and corresponding velocity component for each electron (in PIC and RSM). We find absorption by electron via AHR happens in a very short interval \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Delta \tau$$\end{document}Δτ for less than half a laser period where \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Delta \psi$$\end{document}Δψ remains close to \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\pi$$\end{document}π (necessary condition for maximum laser absorption) and then \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Delta \psi$$\end{document}Δψ drops to its initial \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\pi /2$$\end{document}π/2 (meaning no absorption) after such short-lived AHR. On the contrary, auxiliary magnetic field near the ECR modifies AHR scenario inside the cluster and also helps maintaining the required phase \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Delta \psi \approx \pi$$\end{document}Δψ≈π for the liberated cluster-electron accompanied by frequency matching for ECR/RECR for a prolonged \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Delta \tau$$\end{document}Δτ (which covers 50–60% of the laser pulse through pulse maxima) even after AHR—leading to jump in \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathcal {E}_A}\approx$$\end{document}EA≈ 30–70\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$U_\mathrm {p}$$\end{document}Up. We note that to realize the second stage of enhanced energy coupling via ECR/RECR, the first stage via AHR is necessary.

www.nature.com/scientificreports/ inside the cluster. However, in the few-cycle pulse-regime ∼ 10-fs or below 29,30 , insufficient cluster-expansion forbids LR. But, as the driven electron's excursion amplitude r(t) increases beyond the harmonic regime of the over-dense potential by the laser, its eigen-frequency �[r(t)] drops below ω M . Anharmonic resonance (AHR) occurs in the anharmonic potential when decreasing �[r(t)] of electron meets ω and the electron is promptly ejected out of the cluster with irreversible energy gain [31][32][33][34] . AHR was shown as a strong collisionless mechanism with short pulses by rigid sphere models (RSM), molecular dynamics (MD) and particle-in-cell (PIC) simulations [31][32][33][34][35][36][37] . So called "vacuum heating" 38,39 for LCI is less clear for I 0 <10 18 W/cm 2 . While many experiments demonstrated energetic electrons [12][13][14][15]17,40 with I 0 < 10 18 W/cm 2 for ≈ 780 − 800 nm, and various mechanisms were proposed through analytical models 8,24,[31][32][33][34]37,41 and numerical simulations 8,[42][43][44] to justify experimental findings; still there is no consensus for the maximum energy E A max that an electron can gain (on an average) for a given set of laser and cluster parameters. In fact, historic experiment 13 claiming multi-keV electron energy and its double peak spectrum were later called into question 14,40 ; but values of E A max in those cases were found about ≈ 2.2 times the electron's ponderomotive energy U p = I 0 /4ω 2 (the average energy of a free electron in an oscillating field ∼ √ I 0 sin ωt ). Similarly, various oscillator models 24,[31][32][33][34]37,41 , MD and PIC simulations 33,34,37,41,42,[45][46][47][48][49] showed E A max near 3.2U p . Thus our extensive survey (see Table 1) reveals that in the collisionless regime of LCI, value of E A max mostly remains close to the famous 3.17U p of the laser-atom interaction [50][51][52] ; except in a few cases 43,44,53 where electron's energy around 8U p (or more) were also reported which is imprecise to us. Possibly, collisional events were much active therein. Nevertheless, the primary objective of this work is to increase E A max of cluster-electron far beyond U p . We concentrate in the fascinating 5-fs (fwhm) short-pulse regime of laser ( I 0 > 10 15 W/cm 2 , = 800 nm) interacting with a small cluster ≈ 3 nm where both CA and LR can be ignored and AHR is applicable. By RSM and detailed PIC simulations here, AHR alone is shown to yield E A max 3.2U p similar to earlier works [31][32][33][34]37,41 . We retrieve the phase-difference �ψ between the driving laser electric field and corresponding velocity component for each electron (in PIC and RSM) in the laser polarization; and find that fast generation of electrons via AHR occurs within a tiny interval �τ where �ψ remains close to π (necessary condition for maximum energy absorption rate). This condition �ψ ≈ π holds only for �τ less than half a laser period T = 2π/ω and then �ψ quickly drops to its initial π/2 (meaning no further absorption) after such short-lived AHR. Though remaining laser pulse has adequate supply of energy to the AHR-freed electron, it can't retain finally to conserve the canonical momentum. Therefore, coupling of this unused laser energy to the AHR-freed electron requires a second mechanism which is envisaged here with an ambient magnetic field B ext ; the electron may be energized meeting the electron-cyclotron resonance (ECR) when its cyclotron frequency � c0 = |eB ext /m 0 | = ω . Such ECR requires a facility for strong B ext ≈ 10-20 kT for 800 nm laser. For CO 2 laser (wavelength ≈ 10.6 µm ), the requirement of ambient magnetic field strength for ECR is lowered and B ext ≈ 1-2 kT is sufficient which seems to be feasible. Self-generated (quasi-static) magnetic field of tens of kilo-Tesla are also often measured in the background of laser-plasma experiments and astro-physical conditions. For example self-generated magnetic fields in the range of 20-46 kT have been measured almost two decades ago in laser plasma experiments 54,55 . Magnetic field strengths in the environment of neutron stars and pulsars exceed 10 kT 56 and typically lie in the range of ≈ 10-100,000 kT. Understanding of the origin of energetic electrons in these strong electromagnetic field conditions are of fundamental interest and the present work may explain them. Recent laboratory demonstration of pulsed magnetic fields from sub kilo-Tesla 57,58 and kilo-Tesla to mega-Tesla 59-63 has already renewed interest in laser-plasma 64-66 community and may serve our purpose. From the practical application point of view, energetic electrons produced by LCI via ECR in the presence of an ambient magnetic field can be helpful for the table-top intense radiation sources (such as x-rays) and particle-accelerators which are useful for medical applications. Energetic electrons may help accelerating plasma ions and neutrals 11 via secondary process (e.g., charge exchange, recombination etc.) which are also useful in medical applications and material science.
We show an enhanced E A max ≈ 30-70U p , almost 15-30 fold increase, with an ambient B ext (in crossed orientation) near the ECR even with non-relativistic I 0 ∼ 10 15 − 2 × 10 17 W/cm 2 . Due to relativistic mass increase with increasing energy ( γ m 0 c 2 ), electrons may quickly deviate from the standard (non-relativistic) ECR condition � c0 = ω , but time-dependent relativistic-ECR (RECR) happens with relativistic electron-cyclotron frequency � c = � c0 /γ (t) = ω during the laser pulse driving. The ambient B ext near the ECR not only modifies AHR scenario inside the cluster, it helps maintaining the required phase �ψ ≈ π as well as frequency matching for ECR/RECR for the liberated electron from the cluster in the free space. This phase matching �ψ ≈ π is maintained for a prolonged duration-�τ extending 50-60% of the 5-fs broadband pulse-through pulse maxima even after the AHR, leading to huge absorption E A max ≈ 30-70U p . Here AHR first sets a transverse momentum with which liberated electron is self-injected (no external injection scheme is required) into the remaining laser field where B ext re-orients its momentum and helps energizing it further in the free-space enforcing improved phase-matching �ψ ≈ π and frequency matching for ECR/RECR. However, to realize the ECR/RECR stage (second stage), a transverse momentum of the electron through AHR (first stage) is necessary. Our PIC results are shown to be well-supported by RSM.

Methods
Rigid sphere model of field driven cluster. In the RSM, cluster is assumed as a pre-ionized spherical nano-plasma of radius R and fixed ionic charge density ρ i . RSM has been widely used for LCI [29][30][31][32][33][34]37,67 without B ext . In this work we first include B ext in RSM to study its effects. Ions provide the potential φ(r) with the spacecharge field www.nature.com/scientificreports/ in which electrons interact in addition to the applied laser field ( E l , B l ) and external B ext . Dynamics of an electron obeys where γ = 1/ 1 − v 2 /c 2 = 1 + p 2 /m 2 0 c 2 is the relativistic γ-factor for the electron, m 0 , q, r, v, p are its restmass, charge, position, velocity and linear momentum respectively with m 0 = 1 , q = e = −1 in a.u. Equations (1)-(3) represent a field driven three-dimensional non-linear oscillator. The coulomb part of E sc ∝ r/r 3 restricts its analytical solution, except in some simplified linear case of E sc ∝ r with continuous (plane-wave) laser field only. For example, see direct laser acceleration (DLA) of electrons from an under-dense, pre-formed plasma channel 65,68-72 assisted by auxiliary fields, e.g., magnetic wigglers, static electric and magnetic fields with I 0 > 10 18 W/cm 2 and corresponding normalized vector potential a 0 = √ I 0 /ωc > 1 . To obtain electrons of MeV energies (or higher), the regime of a 0 > 1 is an obvious choice. Such pre-formed plasma channels are very long (typically tens of ) and relativistically intense laser has to propagate several which then sets up electro-static fields in the channel with associated self-generated quasi-static magnetic fields. Electrons are injected into the channel or drawn from the plasma itself and guided by the channel's fields and the applied laser field. If the ambient magnetic field is in the direction transverse to the laser polarization, then energy of electrons can be increased and ECR may happen if such magnetic field satisfies the ECR condition. This work, however, reports other unexplored regime of DLA with I 0 < 10 18 W/cm 2 using short-pulsed light and a constant B ext for an over-dense cluster plasma electrons.
The field E sc imparts oscillatory motion in r , whereas B ext imparts rotation in the plane perpendicular to B ext (in r ⊥ ) to an electron. Combining these two motions, the position dependent squared effective-frequency ω 2 eff [r(t)] of electron in the RSM [using (1)-(2)] can be formally obtained as The term (γ E sc /r + � 2 c0r ⊥ ) ·r ⊥ /γ 2 represents motion due to combined space-charge and v × B ext field in r ⊥ plane and (γ E sc ·r || )/γ 2 r represents motion in r || along B ext . The unit vectors r,r ⊥ indicate frequencies are valid only for motions in those directions. Equation (5) may be regarded as the relativistic extension to its non-relativistic variant 33,34,37,67,73 for c0 = 0 and γ = 1 . When E sc = ω 2 M r and γ ≈ 1 , it gives harmonic oscillator frequency ω 2 eff [r(t)] ≈ (ω 2 M + � 2 c0 ) for low |B ext | values inside the cluster where ω M ≫ � c0 , r || ≈ 0 . It may also be looked upon as upper-hybrid electron frequency 74,75 in magnetized plasmas. Due to non-linear E sc and the relativistic non-linearity imposed by strong E l , B l , B ext ; the ω 2 eff [r(t)] drops from ω 2 eff [0] ≈ (ω 2 M + � 2 c0 ) for increasing r > R . An electron may absorb laser energy by AHR when its ω eff [r(t)] dynamically meets the condition ω eff [r(t)] = ω for increasing r > R above a certain strength of |E l | . Significance of AHR was explained [31][32][33][34][35][36][37]41 .
Here, we shall also present modified AHR (using Eq. (5)) with B ext .
Regeneration of previous RSM results: single electron dynamics with laser field only. We begin with energy absorption and associated electron's dynamical variables as a conventional case of LCI without B ext . Figure 1 shows (column-wise) normalized co-ordinates ( x = x/R, y = y/R, z = z/R ), squared effective frequency ω 2 eff = ω 2 eff /ω 2 , total energy E = ((γ − 1)m 0 c 2 + qφ)/U p in units of U p and corresponding momenta ( p x = p x /c, p y = p y /c, p z = p z /c ) vs time t/T of an initially bound electron [ r(0) = 0, p(0) = 0 ] in the RSM potential, when driven by 5-cycle pulse of I 0 = 7.13 × 10 16 W/cm 2 . Associated fields E x =x.(E l + E sc ), E y =ŷ.E sc , E z =ẑ.E sc and E l =x.E l vs t/T are also plotted, purpose of which will be evident when we consider B ext later. Left panels (a1,b1,c1) and right panels (a2,b2,c2) are the cases with B l = 0 and B l = 0 respectively. www.nature.com/scientificreports/ Without B l (in many works B l was neglected 21,27,31,76 since max |B l | = E 0 /c ≪ 1 ), it is shown that electron starts ( Fig. 1a1) with the binding energy E = −1.5ω 2 M R 2 /U p at t/T = 0 , it oscillates in the potential with increasing amplitude in x (while y = 0, z = 0 ) as the total field E x oscillates in time (while E y = 0, E z = 0 ) and approaches to the peak value E 0 of E l (Fig. 1c1) around t/T ≈ 2 . Inside the potential, for r/R ≤ 1 , E x is suppressed due to opposite phase of x.E l and x.E sc . As long as r/R ≤ 1 , ω 2 eff continues ( Fig. 1a1) at its initial value ω 2 M /ω 2 ≈ 9.1 . When E x increases sufficiently strong (due to reduced phase mismatch between x.E l and x.E sc ) leading to increasing r/R > 1 ; ω 2 eff falls rapidly, it meets the AHR condition ω 2 eff = 1 around t/T ≈ 2.1 (marked by horizontal dashed line and vertical shaded bar) and then electron leaves the cluster forever ( Fig. 1a1) with E > 0 associated with non-zero transverse momentum p x (in Fig. 1b1). After the AHR, E x follows E l . Though LR can not happen, AHR is dynamically met here leading to the electron's removal from the cluster with E > 0 and non-zero p x eventually. Similar results (neglecting B l ) are shown in Refs. [31][32][33][34]37,41 .
Considering B l now, Fig. 1a2,b2,c2 show indistinguishable variation of ( x, y, ω 2 eff , E, p x , p y , E x , E y ) with respective Fig. 1a1,b1,c1. Also z ≈ 0, p z ≈ 0 (in a2, b2) before the occurrence of AHR near t/T ≈ 2.1 , since v × B l field along z is much weaker and leads to a negligible E z [ Fig. 1c2, clearly seen in zoomed inset plots in (b2,c2)]. As the electron is liberated (Fig. 1a2) via AHR around t/T ≈ 2.1 with dominant velocity in x (Fig. 1b2), the v × B l field imparts a forward momentum p z along the laser propagation ( Fig. 1b2) and its z co-ordinate sharply increases (Fig. 1a2) by many times R. Electron is now emitted in the z − x plane with an angle θ ≈ arctan(p x /p z ) in contrast to Fig. 1a1,b1 where electron is emitted only along the polarization axis. Though I 0 < 10 18 W/cm 2 , the liberated electron via the AHR process around t/T ≈ 2.1 is self-injected into the remaining laser pulse with some forward momentum p z > 0 and transverse momentum p x ; and from this time onward electron's acceleration resembles the standard DLA. Clearly, inclusion of B l here yields (Fig. 1a2,b2,c2) different electron dynamics (see also Mulser et al. 31 ) for LCI than neglecting it 21,27,31,41,[76][77][78] in previous works.
However, both the cases in Fig. 1 show maximum attainable energy max E = 8 (marked by upper horizontal dashed line) near the laser peak at t/T = 2.5 ; but the electron retains only a lower value of energy E A = E(τ ) ≈ 2.1 in the end. We may compare these two limits of E with the laser-driven electron-atom www.nature.com/scientificreports/ re-collision model [50][51][52]79,80 of harmonic generation where max E of electron may go up to ≈ 8 during the pulse, but the returned electron when re-collides with the parent ion has a lower E ≈ 3.17 which is often manifested as a harmonic cut-off energy. In laser-cluster experiments, an electron's final energy is reported to be less than the above mentioned laser-atom interaction case [81][82][83]  3.17 in the collision-less case. Thus, though the role of B l can not be neglected for altering the electron dynamics (Fig.1a2,b2,c2) at a I 0 > 7.13 × 10 16 W/cm 2 (where peak magnetic field can be substantial > 2.44 kT), the average E max A 3.17 seems to follow (see Table 1) for the traditional LCI. The aim of the paper is to increase this limit far beyond E max A ∼ 3.17 with an ambient B ext .
PIC simulation. We also study LCI with/without B ext using three-dimensional PIC simulation code 33,34,41,76,78,84,85 . The same deuterium cluster with number of atoms N = 2176 is considered. Atoms are placed in a cubical computational box according to the Wigner-Seitz radius r w ≈ 0.17 nm (giving cluster radius R = r w N 1/3 ≈ 2.05 nm) so that center of the cluster coincides the center of the computational box of side L = 24.6R . Initially laser E l (t) ionizes all neutral atoms D to D + (assuming over-the-barrier ionization, OBI 86 which is valid for I 0 > 10 15 W/cm 2 ) after reaching a critical strength E c = |E l (t)| = I 2 p (Z)/4Z , where I p (Z) is the ionization potential for charge state Z = 1 . Such a fully ionized cluster initially acquires a charge density ρ i /ρ c ≈ 27.87 and ω M /ω ≈ 3.05 at 800 nm. Thus cluster parameters are kept as the RSM. The position and velocity of a newly born electron (after the OBI) are assumed same as the parent atom/ion conserving the momentum and energy. Subsequent movement of more mobile electrons from the relatively less mobile ions by the driving fields create/modify space-charge field E sc (r, t) . Thus E sc (r, t) = −∇φ(r, t) and corresponding potential φ(r, t) in PIC are time-dependent and start from zero contrary to the RSM.
A PIC electron/ion has the same charge to mass ratio of a real electron/ion. The equation of motion of the j|k-th PIC electron/ion (j for electron and k for ion) reads j|k /c 2 , v j|k , r j|k , m j|k , q j|k are relativistic momentum, velocity, position, mass, and charge of a PIC electron/ion respectively. In the present case, values of potential and charge density) with time-dependent monopole boundary condition. Interpolating φ G to the particle position corresponding potential φ(r j|k , t) is obtained. Field E sc (r j|k ) = −∇φ(r j|k ) in (8) is obtained by analytical differentiation 76 of interpolated φ(r j|k ) locally at r j|k . Equation (8) is solved by VVM using laser fields (6)- (7). Total absorbed energy E(t) = l q l φ l + p 2 l /2m l is obtained by summing over all electrons and ions. For the 5-fs pulse (used here) contribution of ion kinetic energy is small and total energy is mainly due to electrons. The numerical parameters in the PIC simulation (spatial and temporal resolution, grid size, number of PIC particles/cell etc.) are carefully chosen for negligible artificial numerical heating. Typically, we have chosen 64 × 64 × 64 grid points (cells) with uniform grid size x = y = z = 16 a.u., time step t = 0.1 a.u., and approximately 15 particles/cell. Two important upgrades are made in the current PIC version: relativistic particle mover based on (i) Runge-Kutta 4-th order method (RK4) and (ii) VVM. It is found that VVM leads to better energy conservation and less numerical heating even for a bigger t than RK4, particularly for the relativistically intense driving fields. Electron-ion collisions are neglected in the current work due to high field values.

Results
New RSM results: electron dynamics with laser and auxiliary B ext . Results in Fig. 1 (right column) show that for the chosen E l , B l configuration, a liberated electron from cluster may also gain a mild forward momentum p z after AHR. The energy-momentum relation p z − p z0 = (γ − γ 0 )c for DLA (without space-charge) suggests that to improve energy gain by the electron, its p z should be increased from the initial p z0 = γ 0 c . Though magnetic field does not work, an auxiliary B ext helps bending electron's trajectory. It may also improve p z of the freed electron. Figure 2a1,b1,c1 show results with B ext = |ẑB ext | ≈ 6.68 kT (corresponding � c0 = ω/2 ) along z for the same (E l , B l ) as in Fig. 1a2,b2,c2 which is considered as a reference. Noticeably, variation of dynamical variables are now very different from the corresponding Fig. 1a2,b2,c2; but the final retained energy of the electron is still E A ≈ 2.1 . The ω 2 eff starts at ≈ (ω 2 M + � 2 c0 )/ω 2 , monotonically drops and passes the AHR line ω 2 eff = 1 (horizontal dashed line) at a little early time t/T ≈ 1.85 (vertical shaded bar) following Eq. (5) contrary to its short-time oscillatory nature just before the occurrence of AHR (Fig. 1a2) near t/T ≈ 2.1 . The vanishing of oscillatory nature of ω 2 eff (Fig. 2a1) and its smooth passage through the ω 2 eff = 1 line is due to additional induced fields E y , E z [though still weak, Fig. 2c1] due to strong B ext =ẑB ext leading to swirling motion in (x, y) inside the cluster similar to the driving by a circularly polarized laser field 33,73 . Thus an external B ext =ẑB ext may modify electron dynamics inside the cluster and the AHR scenario. The ( x, y, p x , p y ) dynamics of the liberated electron tends to follow cyclotron motion; both p z and E grow up-to a maximum (note that max E ≈ 10.5 now exceeds the conventional E = 8 line without B ext ) near the pulse peak. But p z drops later (Fig. 2b1) leading to lesser (Fig. 2a1) final energy E A = E(τ ) ≈ 2.1 as in Fig. 1a2 though electron dynamics drastically differ from Fig. 1a2,b2,c2.
With a higher B ext = |ẑB ext | ≈ 13.37 kT corresponding to � c0 = ω (ECR), Fig. 2a2,b2,c2 show a significant jump in the final absorbed energy upto E A ≈ 36 (far exceeding the conventional E max A ≈ 3.17 ) associated with a jump in the corresponding final p z ≈ 0.3 . Most of the arguments relevant to Fig. 2a1,b1,c1 apply here also. Additional inset plots are zoomed view of dynamical variables near AHR and inside the cluster. Due to higher www.nature.com/scientificreports/ B ext =ẑB ext , induced fields E y , E z in the cluster (Fig. 2c2) are also marginally stronger, AHR scenario is marginally modified and note that, even in this case p z , z are almost zero inside the cluster. Distinctly, after the AHR near t/T ≈ 1.85 , liberated electron follows almost exact cyclotron motion in the x − y plane (evident from x, y, p x , p y variation) due to the stronger B ext =ẑB ext , while its p z and E continuously increase to p z ≈ 0.3 and E ≈ 36 during t/T ≈ 2 − 4 followed by saturation, though laser field envelope (Fig. 2c2) weakens after its peak.
Thus an external magnetic field -assisted electron acceleration from a laser-driven cluster is shown to enhance electron's energy by 10-12 times than the conventional limit of E max A ≈ 3.17 , particularly near the ECR frequency � c0 = ω . This encouraging new result needs further investigation.
Above results with B ext show laser absorption happens mainly in two stages. In the first stage, electron undergoes AHR (may be modified by B ext ) and comes out of the cluster with low positive energy and non-zero transverse momentum. Later, in the second stage, it is fully controlled by the remaining E l , B l and B ext with an increase in E , i.e., absorption rate (see Fig. 2a2) for t/T ≈ 2 − 4 . This second stage may be termed as magnetic field assisted DLA. However, to realize the second stage, energy absorption by electron and its liberation from the cluster in the first stage is necessary, otherwise v × B ext fails.
The temporal phase dynamics. Equation (4) implies that rate of absorption dγ m 0 c 2 /dt by an electron approaches to zero (or negligible) for phase angle �ψ between its velocity and the corresponding driving electric field approaching to the odd integral multiple of π/2 . From Eq. (4), one may apparently conclude no role of B ext for enhanced absorption in Fig. 2a2,b2,c2. Note that in the second stage of energy absorption, where role of E sc is nil, Eq. (4) simplifies to dγ m 0 c 2 /dt = qv x E x = qv x E l ; and, though B ext can not alter E l , it may re-orient the phase of v x (see the cyclotron orbit) w.r.t. E x . To probe this underlying physics, we numerically retrieve phase angles ψ v x , ψ E x , ψ E l w.r.t. central frequency ω of respective v x , E x , E l (since components along the laser polarization matter the most) for two cases: (i) with E l , B l only for Fig. 1a2,b2,c2 and (ii) with E l , B l and B ext ≈ 13.37 kT for Fig. 2a2,b2,c2. Those ψ v x , ψ E x , ψ E l and �ψ = |ψ v x − ψ E x | vs time are plotted in Fig. 3 (see caption). Vertical shaded bars are the respective AHR regions (see Figs. 1, 2) after which the electron is mostly free from spacecharge fields of the cluster and respective ψ E x goes hand in hand with ψ E l in both cases. Little deviation of ψ E l from π/2 far away from the pulse center (at t/T = 2.5 ) is due to 5-cycle broad-band pulse (ideally it is π/2 for  Fig. 1. Panels in the left (a1,b1,c1) and right (a2,b2,c2) columns are with B ext ≈ 0.028, 0.0569 a.u. corresponding to non-resonant � c0 /ω = 0.5 and resonant � c0 /ω = 1 (ECR case) respectively. In the ECR case AHR occurs little early around t/T ≈ 1.85 , and E A reaches up to E(τ ) ≈ 36 compared to E A = E(τ ) ≈ 2.1 in Fig. 1a2; corresponding momenta and excursion also significantly vary after the electron is freed via AHR near t/T ≈ 1.85 . Inset plots (in a2,b2,c2) show zoomed view of dynamical variables near AHR and inside the cluster due to strong B ext . Other laser and cluster parameters are as in Fig. 1 www.nature.com/scientificreports/ a monochromatic pulse ∼ sin ωt ). Respective ψ v x , ψ E x in (i) do not differ from (ii) and �ψ ≈ 0.5π remains upto t/T ≈ 1.4 . After this time, 9π is maintained for a tiny duration �τ ] followed by its gradual drop through the AHR region and saturation near �ψ ≈ 0.5π afterwards. On the contrary, in (ii) an instantaneous phase swing occurs (near t/T ≈ 1.4 ) for ψ E x = −π → π after quick dropping to −π . Later, though ψ E x → ψ E l ≈ 0.5π , the phase ψ v x is dynamically tilted in a way that a value of �ψ ≈ π is brought about by the auxiliary B ext for a long duration t/T ≈ 1.75 − 3.0 (leading to high absorption rate in Fig. 2a2) from pre-AHR to post-AHR time through the pulse maxima; then �ψ gradually drops as ≈ π → 0.5π for t/T ≈ 3.0 − 4.0 where absorption slows down and finally saturates at a higher E A ≈ 36 in Fig. 2a2. Thus an auxiliary B ext near the ECR helps maintaining the required �ψ ≈ π for enhanced laser absorption in the second stage.

Non-interacting multi-electrons in RSM.
A single-electron dynamics (as studied by RSM above) is important to understand the physics of LCI, but can not answer some other aspects, e.g., fraction of electrons leaving the cluster (outer ionization fraction) and their energy distribution. In the single-electron case, outer-ionization fraction assumes only 0,1 (electron is either inside or outside the cluster). In a real system, however, some electrons may remain bound and outer-ionization fraction may attain any value between (0,1) depending upon laser and cluster parameters. A single-electron case may over-estimate/under-estimate electron energy compared to the realistic multi-electron case where per-electron energy may be averaged out. Moreover, different electrons become free from the cluster at different times, and participate in the magnetic field assisted DLA differently. To answer these aspects we distribute all N = 2176 electrons inside the cluster randomly (or uniformly) to mimic a multi-electron system by RSM where electrons are assumed non-interacting among them. For brevity, we compare these multi-electron results of RSM along with detailed PIC simulation in the following section where particle-particle interactions are taken care self-consistently.
Absorption studies with PIC simulation and comparison with RSM. Figure 5 compares time vs average energy (per electron) between PIC and RSM results for B ext = 0 and |ẑB ext | = ω (ECR case) at  Fig. 1a2,b2,c2 without B ext and (ii) in Fig. 2a2,b2,c2 with B ext = |ẑB ext | ≈ 0.0569 (ECR case). Phases are numerically calculated by FFT w.r.t. the central frequency ω with a sliding time-window duration T = 2π/ω . In (i) relative phase �ψ = |ψ v x − ψ E x | , on an average, stays close to π (or 0.75π ) for a short-while �τ < T/2 . In (ii) �ψ stays close to π for a prolonged duration �τ > T through the pulse peak (before falls to π/2 at t/T ≈ 4 ) leading to E A ≈ 36 in Fig. 2a2.  Fig. 4 and non-interacting multielectrons (RSM-MP) as described in section above are also included. RSM-SP over-estimates the RSM-MP case for final energy E A = E(τ ) when B ext = 0 , but PIC result ( E A ≈ 0.5 ) follows RSM-MP more closely. For the ECR case, however, E A ≈ 36 in PIC remains little higher than RSM-MP, which is due to early ejection of electrons with non-zero transverse momentum via AHR from the self-consistently developing potential and electrostatic restoring fields (starting from zero) in PIC. Note that E(t) starts increasing one-period earlier ( t/T ≈ 1.2 ) in PIC than the RSM and so as the ECR for those early leaving PIC electrons. Almost 60 − 70 fold increase in E A ≈ 0.5 → 36 is obtained in PIC and RSM-MP due to |ẑB ext | near ECR. S c a n n i n g t h r o u g h r a n g e o f v a l u e s o f |ẑB ext | ≈ (0 − 2ω) , f o r d i f f e r e n t I 0 ≈ 1.83 × 10 15 W/cm 2 − 1.83 × 10 17 W/cm 2 and same 5-fs pulse duration, results in Fig. 6 are obtained by PIC and RSM-MP in the end of the pulses. At a low intensity (Fig. 6a) absorption peak (at 26 U p ) occurs almost at the ECR condition � c0 = ω (vertical dashed line) as clearly exhibited by PIC simulation where electrons can undergo AHR at ease, become free with a transverse momentum for the ECR in the next stage; whereas RSM-MP shows almost zero absorption since AHR is not met (first stage fails) and electrons can't be freed from RSM potential with a transverse momentum at this low intensity (RSM greatly under-estimates absorption here,  Comparison of PIC and RSM results: average total absorbed energy E(t) = E(t)/NU p per particle in U p vs t/T with B ext = 0 (dashed lines, conventional case of Fig. 1a2,b2,c2) and with B ext = ω (solid lines, ECR case of Fig. 2a2,b2,c2). RSM results with single-electron (RSM-SP) and non-interacting multi-electrons (RSM-MP) both justify PIC results. Though absorption starts early in PIC, final absorbed energies E A = E(τ ) with/without B ext are comparable with the RSM cases. www.nature.com/scientificreports/ since its needs a threshold intensity 31,32 ). Therefore, as the peak intensity increases, absorption peaks show-up gradually (Fig. 6b,c) for RSM-MP due to gradual removal of electrons via AHR (preferably) from surface to the cluster center, but ECR absorption peak occurs always for PIC. Finally, at a higher I 0 ≈ 1.83 × 10 17 W/cm 2 , PIC and RSM (almost overlap) show very good quantitative agreement in Fig. 6c. Absorption peaks ≈ 65U p , 45U p in Fig. 6b,c give average energy E A ≈ 0.27, 0.49 MeV. The gradual right-shift of the absorption peak from the ECR condition � c0 = ω (vertical dashed line) with increasing intensity is due to the relativistic modification of � c = � c0 /γ in dipole-approximation. Since γ is time-varying (during the pulse) and different for different electrons, the time-dependent relativistic-ECR occurs for electrons when � c (t) = � c0 /γ (t) = ω (call it RECR, see Fig. 7). It emphasizes quick slippage of electron from the RECR condition as soon as its γ (t) > 1 . Therefore, to satisfy the RECR for γ > 1 , a higher c0 (or higher B ext ) is required-as manifested by gradual right-shift of the absorption peak (Fig. 6-c) with increasing intensity. Moreover, laser pulse being broadband with frequencies ω, (1 ± 1/n)ω , RECR may happen in a wider frequency range and contribute to broadening of resonance-width about the absorption peak in Fig. 6. Energy is shown normalized by corresponding U p (left y-axis) and in atomic units (right y-axis). At a low intensity (a) absorption peaks almost at the ECR condition � c0 = ω (vertical dashed line) as clearly exhibited by PIC where electrons undergo AHR at ease and become free with transverse momentum for the ECR in the next stage; whereas RSM-MP shows almost zero absorption since AHR is not met (first stage fails) in RSM. As I 0 increases to moderate values (in b,c) absorption peaks show-up in RSM-MP due to meeting of AHR followed by ECR. For high intensity RSM-MP justify PIC results quantitatively. Gradual right-shift of the absorption peak from ECR condition � c0 = ω with increasing I 0 is due to relativistic modification of � c = � c0 /γ for γ > www.nature.com/scientificreports/ Frequency and phase dynamics of PIC electrons. To elucidate further, we retrieve the relativistic anharmonic eigen-frequency eff and cyclotron-frequency c for each k − th PIC electron as (see also Eq. (5)) Figure 7 shows temporal variation of � eff /ω (green, left y-axis) and � c /ω (red, right y-axis) of cluster electrons for |ẑB ext | ≈ 0.02, 0.0569, 0.07 a.u. at I 0 = 7.13 × 10 16 W/cm 2 (left column, a1,b1,c1) and |ẑB ext | ≈ 0.02, 0.0569, 0.078 a.u. at I 0 = 1.83 × 10 17 W/cm 2 (right column, a2,b2,c2) corresponding to PIC results in Fig. 6b,c respectively. Chosen values of |ẑB ext | for each intensity represent data points A, B, C (at the (9) � 2 eff =r · E sc (r k )/γ k r k  Fig. 6b,c respectively. Vertical shaded region indicates AHR region where � eff /ω of each electron starts from zero, reaches different maximum, then drops to zero passing through AHR when electron is freed from the cluster potential with excursion r/R ≫ 1 (gray) and non-zero transverse momentum. Horizontal dashed lines represent frequencies of the broadband pulse where ECR/RECR are expected. At low B ext values, ECR is not met (a1,a2), laser absorption is mainly due to AHR occurring for t/T ≤ 2 (1st stage, vertical shaded region). As B ext increases, γ of electrons increase, all frequencies of the broadband pulse gradually come under ECR/RECR condition (second stage) with decreasing � c /ω as one passes (b1,b2) to (c1,c2). In (c1,c2) ECR/RECR is hit around the peak of the pulse (at t/T = 2.5 ) with central frequency ω as well as with side-bands 1.2ω, 0.8ω leading to higher absorption in (c1,c2) compared to the case (b1,b2). In (b1,b2) ECR is hit in the beginning of the pulse with ω when laser field is relatively weak, then RECR with the sideband at 0.8ω near the pulse peak and in the pulse end (for b2). Note that ECR/RECR occurring at very early time ( t/T < 1.5 ) or very late time ( t/T > 4 ) are less effective due to weak laser field. Other laser and cluster parameters are as in Fig. 1. See also Fig. 8  www.nature.com/scientificreports/ tail, at the non-relativistic ECR condition � c0 = ω , and at the peak) on the PIC absorption curves in Fig. 6b,c. Normalized position r/R of electrons (gray) show their distances w.r.t. center of the cluster. Corresponding phaseangles ψ v x , ψ E x of respective v x , E x along the laser polarization and the phase difference �ψ = |ψ v x − ψ E x | for each k − th PIC electron are computed by FFT (as in RSM Fig. 3) and shown in Fig. 8 for both intensities. Wavy dashed lines (Fig. 8) indicate respective phase angles for average values v x /N, E x /N of all N electrons showing average system behavior. Contrary to the RSM, � eff /ω of each PIC electron starts from zero and reaches different maximum (Fig. 7) when its r/R drops towards the potential minimum. Then � eff /ω each PIC electron drops to zero passing through AHR 33,34,37 similar to the RSM and electron is freed from the cluster potential with r/R ≫ 1 . Shaded vertical bar (in Fig. 7) highlights this AHR dominated region (1st stage) during initial time of the laser pulse. Since different electrons undergo AHR at different times and comes out with different non-zero transverse momentum, the exact extent of the 1st stage and the beginning of 2nd stage (ECR stage) with B ext is difficult to draw (i.e., minor overlap happens and 2nd stage starts early for early leaving electrons via AHR) with all electrons together. However, from the vanishing of � eff /ω → 0 and increasing r/R ≫ 1 it is clear that AHR domain (1st stage) is mostly limited below t/T ≈ 2 for I 0 = 7.13 × 10 16 W/cm 2 ( t/T ≈ 1.6 for I 0 = 1.83 × 10 17 W/cm 2 ) and shrinks with increasing intensity. At low B ext values, ECR is not met, �ψ quickly falls to π/2 after initial rise towards π due to AHR (mostly occurring) for t/T ≤ 2 . As B ext increases towards ECR, �ψ is gradually lifted towards π , and it is maintained for a longer duration �τ ≈ 60 − 70% of the pulse through pulse maximum leading to higher absorption in Fig. 6b,c even after AHR. Phase angles are numerically computed by FFT as in RSM Fig. 3. See also Fig. 7 for corresponding frequency dynamics. Other laser and cluster parameters are as in Fig. 1 www.nature.com/scientificreports/ At low |ẑB ext | ≈ 0.02 (or without it) as in Fig. 7a1,a2, the frequency matching for ECR can not happen, the phase difference �ψ continues to π/2 after initial rise towards π shown in respective Fig. 8a1,a2 due to shortlived AHR occurring below t/T < 2 . Hence absorbed energy remains low ( < 3U p ) without initiating the second stage. In these cases not all electrons are freed (Fig. 7a1,a2), many of them may comeback inside the cluster later time, and may be liberated again through another AHR, e.g., see after t/T > 4.

Scientific Reports
As B ext increases (see caption of Fig. 7), γ of electrons increase, all frequencies of the broadband pulse (shown by horizontal dashed lines) gradually come under ECR/RECR condition with decreasing � c /ω as one passes Fig. 7b1,b2-c1,c2; accompanied by gradual lifting of �ψ towards π even after AHR with more time elapsed near π as in respective Fig. 8b1,b2-c1,c2. Also, as RECR is met with the central frequency ω near the pulse peak (Fig. 7) at t/T = 2.5 and respective �ψ is maintained near π for a longer duration �τ ≈ 60 − 70% of the pulse through pulse maximum (Fig. 8), it leads to higher absorption ( E A ) in Fig. 6c, b. Thus, not only frequency matching � c /ω = 1 for ECR/RECR is satisfied, the required phase matching condition �ψ ≈ π is also simultaneously satisfied by PIC electrons for all cases in Fig. 6 (same are checked with electrons in RSM-MP for Fig. 6, not repeated) for enhanced absorption peak about 30-70U p .

Discussion and summary
We study laser-deuterium cluster interaction with short 5-fs (fwhm) laser pulses ( I 0 > 10 15 W/cm 2 , = 800 nm) in presence of external magnetic field B ext ≈ 10-20 kT using RSM and three-dimensional PIC simulations. For the standard case, without B ext , our extensive survey on laser-cluster interaction finds that average energy per electron E A most often remains around 3.2U p or less. Without B ext , first we show that AHR alone may yield E A 3.2U p similar to earlier works 24,[31][32][33][34]37,41,42,[45][46][47][48][49]67 even with the inclusion of the laser magnetic field B l . We then retrieve the phase-difference �ψ between the driving laser electric field and corresponding velocity component for each electron (in PIC and RSM) in the laser polarization and find that generation of electrons via AHR occurs within a short interval �τ where �ψ remains close to π (necessary condition for maximum energy absorption rate) only for �τ less than half a laser period. After that �ψ quickly drops to its initial π/2 , leading no further absorption. Thus AHR is found to be very short-lived. Though remaining laser pulse supplies energy temporarily, the AHR-freed electron can not retain finally. Therefore, coupling of this unused laser energy to the AHR-freed electron is envisaged here through a second mechanism with B ext , namely ECR, when electroncyclotron frequency � c0 = |eB ext /m 0 | = ω.
We show an enhanced average energy per electron E A ≈ 30-70U p with an ambient B ext (in crossed orientation) near the ECR even with non-relativistic I 0 ∼ 10 15 − 2 × 10 17 W/cm 2 . Due to relativistic mass increase with increasing kinetic energy ( γ m 0 c 2 ), electrons quickly deviate from the standard (non-relativistic) ECR condition � c0 = ω , but time-dependent relativistic-ECR (RECR) happens with relativistic electron-cyclotron frequency � c = � c0 /γ (t) = ω during the laser pulse driving. The ambient B ext near the ECR not only modifies AHR scenario inside the cluster, it also helps maintaining the required phase �ψ ≈ π as well as frequency matching for ECR/RECR for the liberated electron from the cluster in the free space for a prolonged duration �τ . We find that �τ extends ≈ 50-60% of the 5-fs broadband pulse -through pulse maxima even after the AHR-leading to huge absorption E A ≈ 30-70U p . Here AHR first sets a transverse momentum with which liberated electron Table 1. Approximate value of maximum average absorbed energy E A max (in units of ponderomotive energy U p ) of an electron from traditional laser-cluster interaction in various published works. These are calculated from the available data/graphs. In most cases E A max remains close to 3.2U p or below.

Data availability
The data that support the plots and findings of this paper are available from the corresponding author on reasonable request. However, due to other novel findings, authors won't be able to make the raw data public.